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1. Introduction 



In this paper we focus on an inverse eigenvalue problem for unitary llessenberg matrices with pos- 
itive subdiagonal elements. Throughout this paper all llessenberg matrices are upper llessenberg 
matrices. This class of matrices bears many similarities with the class of Jacobi matrices, i.e. real 
symmetric tridiagonal matrices with positive subdiagonal elements. Any matrix in either class is 
normal and has distinct eigenvalues. Both n x n Jacobi matrices and n X n unitary llessenberg 
matrices with positive subdiagonal elements can be parameterized by 2n — 1 real parameters. This 
is obvious for Jacobi matrices; for unitary llessenberg matrices this parameterization is described 
below. Since unitary llessenberg matrices with positive subdiagonal elements are determined by 
0{n) parameters, one can develop efficient algorithms for this class of matrices. These algorithms 
are analogous with algorithms for Jacobi matrices. For example, the unitary QR algorithm, de- 
scribed in [Gr2], has many similarities with the QR algorithm for Jacobi matrices. Another example 
is provided by the divide-and-conquer methods that have been developed for both the tridiagonal 
and unitary eigenproblems [Cu], [DS], [GRl], [GR2], 

In this paper we show another analogy of unitary llessenberg matrices with Jacobi matrices, namely, 
that a unitary llessenberg matrix H with positive subdiagonal elements is uniquely determined 
by its eigenvalues and the eigenvalues of a unitary rank-one perturbation of H . The matrix H 
can be constructed using 0{iv) arithmetic operations using an inverse unitary QR algorithm. 
Similar results for Jacobi matrices are well-established [BGl], [BG2], [GH]. The inverse unitary QR 
algorithm is analogous with the algorithm described in [GH] for Jacobi matrices. 

2. Unitary Hessenberg Matrices and Szego Polynomials 

We refer to a finite Schur parameter sequence of length n as a sequence of complex numbers {7j}j=i 
with | 7 j| < 1 for I < j < n and | 7 ^| = 1. Also define the complementary Schur parameters 
by Gj \/l — |7j|“. Associated with the finite Schur parameter sequence is a unitary 

Hessenberg matrix H with positive subdiagonal elements 

H = /f(7i, . . . ,7n~l?7n) ♦= ( 71 )^ 2 ( 72 ) . . . G'n_i(7n_i)(5n(7n), (2.1) 

where the Givens reflector Gj{^j) is the identity matrix of appropriate size except for the 2x2 
principal submatrix 



j i + 1 




1 

1 

1 


j i + 1 







with the bar denoting complex conjugation. The matrix on the right in the product (2.1) is defined 
by Gn(^n) •= diag[l, 1, . . . , 1, — 7 n]. The nonzero entries of H = [Vj,k]]^k=i given by 

Vj+iJ *= and rjj^k := for 1 < j < k, where 70 := 1. 

It is easy to see that every n X n unitary Hessenberg matrix H = with > 0 is 

uniquely determined by a finite Schur parameter sequence of length n. In fact, the Schur parameters 
{ 7 j}j=i and the complementary Schur parameters ^an be determined from H by 

l<j<n; 7j = -77ij/(7ia2 • • 1 < i < n. 

Hence, we have a one-to-one correspondence between Schur parameter sequences of length n and 
nXn unitary Hessenberg matrices with positive subdiagonal elements. This Schur parameterization 
of unitary Hessenberg matrices with positive subdiagonal elements shows that these matrices are 
determined by 2n — 1 real parameters: the real and imaginary parts of 7 j for 1 < i < n, and the 
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argument of 7n- To avoid numerical instability, however, we also retain the coinj)leinentary Schur 
parameters. 

Let II = Ilfi •= //('/i , . . . , 7 n-i ^ 7 n)^ II k := (71 )Go(72 ) • • • )6\.(7;J be the 

leading principal submatrix of II = //„ of order k. Introduce the functions 

<^,.(A):=ef(A/-//,)-'ei, 
rl-k(X) := dci{XI - Ih), 

7t,(A) := det(A/-7/"_i). 

where Cj := [1, 0, . . . , O]^, and denotes the trailing principal submatrix of 7/ of order — 1. 

Thus, <^a-(A) = 7T;t(A)/t/7-(A). The following proposition can be verified by induction. 

Proposition 2.1. The polynomials V^7-(A) and 7T;^(A), k > 0, satisfy 



l/’<:(A) TTkiX) 




X Ik 


^ 

1 

1 

1 




V’o(A) 7To(A) 




1 1/A 


^'k{X) T^k{X) 




ikX 1 






v'-o(A) ^-o(A) 




1 0 



It follows that for each k\ the i)olynomials 4^k{^) A^^t/7.(1/A) are obtained by reversing and 

conjugating the coefficients of V’/:(A). From tlie initial conditions we recognize V7-(A) to be the kxh 
Szego polynomial determined by the Schur parameter sequence {7y }^_j. □ 

The Szego polynomials are orthogonal with respect to a discrete measure on the unit 

circle. This measure assigns a positive weight to each zero of fn(A). These weights are the 
numerators in the partial fraction decomposition 



</>n(A) 



n 



E 



A- A,/ 



see [Grl] for details. 

Let II = with A = diag[Ai , A2, . . . , A^^] and U unitary, be the spectral resolution of 

II = //(71,.. .,7n-i57n)? where + denotes transposition and complex conjugation. Let u := r^ci 
be the vector containing the first components of the eigenvectors of 77 . Since 77 has nonzero sub- 
diagonal elements, every entry of u = [ui, V2 , .... Vn]^ is nonzero, and we normalize U so that each 
> 0. Then 

n o 

<;^>n(A) = U*(A 7 - y . . 

k=\ ^ 

Hence, the weights u-7- = vl are guaranteed to be positive, and ^'k = F 



3 . The Inverse Unitary QR Algorithm 

Given n distinct unijnodular complex numbers {A;:}5!_j and associated positive weights 
we can construct a unitary Hessenberg matrix II with the nnd 17- equal to the eigenvalues and 
first components of the corresponding eigenvectors of 77 , respectively. This construction is achieved 
using an inverse QR algorithm, which is analogous with the ])rocedure of [GH] for real symmetric 
tridiagonal matrices. 



The required Hessenberg matrix is obtained by performing a seciuence of elementary unitary simi- 



laritv transformations to transform the matrix 



u 

A 



to a Hessenberg matrix 



I 



II 



without 
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using or changing the arbitrary entry 8, Then II = //( 71 , . . .,7n~i,7n) has the desired eigenvalues 
and associated eigenvectors. 

The idea is to build the Hessenberg matrix by adding weight-abscissa pairs one at a time. Suppose 
that we have constructed the unitary Hessenberg matrix Hm •= //( 71 , . . . , 7 ^- 1 , 7m), for some 
m < n, corresponding with the weight-abscissa pairs Let ao :=: and 

assume that the first components of the eigenvectors of Hm axe {(vk/(^o)}^-i» order to add the 
weight-abscissa pair A) and construct the corresponding (m + 1) x (m+ 1) unitary Hessenberg 
matrix := H{X[^ Ao, . . . , we perform a sequence of unitary similarity transformations 

to put the {m + 2) X (m + 2) matrix 






6 

V 

(^oei 



V aoe^ 

A 0" 

0 Hm 



S V ao 
V A 

ao X X X X X 

X X X X X 

X X X X 

XXX 
X X- 



into Hessenberg form without changing 6, Let Gq := and Oq := — v/ag. Then 



G2(oo)/7^^^ =: 



■ <5 V ao 

ag X X X X X 

0 X X X X X 

X X X X 

XXX 
X X 
X 



X 
X 
X 
X 
X 
X - 



is a Hessenberg matrix with a trailing principal (m+ 1) x {m + 1) submatrix, which is both unitary 
and of Hessenberg form. On the completion of the similarity transformation of H^^\ we obtain 



H^-'>G;{ao) 




X X X X X 

X X X X X 

® X X X X 

XXX 
X X 
X 



X 

X 

X 

X 

X 

X - 



(3.1) 



The circled element in (3.1) forms a “bulge”, wdiich is to be chased down along the subdiagonal in 
order to obtain a matrix of Hessenberg form. Define ^ 3 ( 01 ) so that G' 3 (ai)//^“^G 2 (ao) =: H^^^ is 
a Hessenberg matrix. Then //(^^G 3 (ai) has a bulge, which we annihilate by multiplying from the 
left with G 4 (q 2 ). Proceeding in this manner, we ultimately chase the bulge off the bottom of the 
matrix, and obtain the Hessenberg matrix which is unitarily similar to H^^\ 

The traihng principal (m + 1) X {m + 1) submatrix of is unitarily similar to 

a unitary Hessenberg matrix with positive subdiagonal elements. The latter matrix is the desired 
Hessenberg matrix H^_^^ = //(7'1, 72^ • ♦ • , 7m+i)* 

This procedure for adding a w^eight-abscissa pair to //^, if implemented by directly manipulating the 
elements of the matrices for 1 < k < 777 , would require 0(7?7^) arithmetic operations. How^ever, 
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wo note that for each k the trailing (?n + l )x(7/i + l ) principal snbinatrix of jg unitary and of Hes- 
senberg form, and, therefore, is unitarily similar to a unitary Hessenberg matrix, denoted by 
with positive subdiagonal elements. Hence, we can carry out the similarity transformations by ma- 
nipulating the Schur parameters of the matrices + 1 Im)- 

This gives rise to a method that requires only 0 ( 7 ?i) arithmetic operations in order to add a weight- 
abscissa pair to Urn- ^Ve refer to this method as the inverse unitary Qli algorithm because of its 
relationship with the unitary QR algorithm presented in [Gr 2 ]. 

Inverse Unitary QR Algorithm: adding a weight-abscissa }>air (r',A), where |A| = 1 . r > 0 . 
CTo := \/crl + v-\ 

ao := -v/oq; /3q := 

for k := 1 , 2, . . . , 777 

•— + |o p; 

^k ’»= Pk-l^{^^k-\ -\-lk^^ '^k-\)/^k'> 

(5k := (5k-\(^klc^k'^ 

7m-|-l 

□ 

Thus, the inverse unitary QR algorithm can be used to construct the unitary n X n Hessenberg 
matrix II = //(71 , . . . , 7 „_i , 7 „) from its eigenvalues and the first components of its normalized 
eigenvectors in 0 ( 77 ") arithmetic operations. 



4, An Inverse Spectral Problem 



Let II — Lf(7i?*‘‘>7n — i’7n) and (j5yj(A) — 



and consider the polynomials 



^n(A) _ L^k 

V’n(A) " A- A,' 



Let a := for some 0 < r < 2r, 



\;,-(A) := (1 - a)A7T;,(A) + qva-(A), 0 < A’ < n. 



(Tl) 



Proposition 4.1. The polynomials {\/:(A)}5!_o are the iiionic Szego polynomials corresponding 
with the Schur sequence {(^lk}k=\- 

Proof, The definition (4.1) and Proposition 2.1 yield 



'uo^y 




A 


Ik 


■\fc-i(A)' 




■\o(A)' 




’ 1 ' 








1 _ 


.a-i(A). 




.s^o(A). 




a 



where ^k{^) •= (1 — <a)A7r/,(A) + ai/^-f A). On setting \^-(A) = a^^-(A), we obtain 



■\a(A)‘ 




A 


Cilk 


■\^-i(A)' 




■\o(A)' 




T 


.U(A)^ 




aifcA 


1 


.U-i(A). 




,\o(A)_ 




1 



An immediate consequence of Projiosition 4.1 is that each zero of \n(A) has unit modulus. W hen 
a = -1 the i)olynomials {\k{^))k=o known as the Szego polynomials of the second kind 
corresponding with {lk](:=i- 



Proposition 4.2. The zeros {/U'}^=o Xn(A) strictly interlace the zeros of 'ipnW on the 

unit circle. Moreover, riLi(/'^7^fc) = «• 

Proof, Let A := for some 0 < 0 < 27T, and let a =: for some 0 < r < 27 t. Then 



^n(A) := 



An(A) _ A — gA^. 

^n(A) A - At 



= 



k=l 



1 — AAt + q( 1 — AAt) 
|A-At|2 



Re(e-/2(l-AAt)) 



= y^2e-'-wt 



k=l 



|A-At|2 



Thus, Im(e = 0 for 0 < 0 < 27 t. Let A^- =: . We may assume that 0 < 6i < Oo < 

• • • < On < 27t. Then 






cos(r/2) — cos(r/2 + Ok — 0) 
"(1 - cos(6>t - 6>))- + sin*(0fc - 0) 



cos(r/ 2 ) — cos(r/2) cos(^t — ^) + sin(r/2) sin(0t — 0) 
= 



1 — cos(0k — 0) 

- (cos(r/2) + sin(r/2)cot((6»t - O)/!)) . 

k=l 

From 0 < r < 27 t and > 0 for all A:, it follows that 

^ sin(r/2) ^wt/sin^((0t - 0)/2) > 0 for 0 # ^t, 1 < k < n. 



k=l 

n 



k=l 



Thus, ”OC as 0 \ 0/^^ and —^ooasOy This shows that 0 — ^ has 

precisely one zero in Consequently, the zeros f.ij of $n(A) strictly interlace the Aj on the 

unit circle. The second statement of the Proposition follows from the fact that yn(d)/'0n(d) = o. □ 

Let Xk and =: for 0 < 0/:,/^/: < 27 t. Then we have ~ 0^) = r. We may 

assume that the arguments have been ordered so that 

0 < 01 < />! < 00 < • • • < l^n-l < On < l^n < + 27T. (4.2) 

Proposition 4.3. With the above notation, and the ordering (4.2), the weights ujk are given by 



^k = -T 



n sin((i/j - 0k)l2) 
_1 jfi 

i=i 



Proof. We have 



. ^ Jli^k-kj) 

/ 1 \ \ T ^ \ \\A^i\^/ \n\^k) j — 1 

(1 - a)A,.. = ^ln,,_(A - A.)^ = ^ 

i=i 
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since \n('^) monic. Tluis. 



n(r‘^*-e -0 n 5>i--'((''j-^0/2) 

,, _ j~j _ _ _OyV*'^/-izi 

[l CX)^^ € ^ ^ 



(^1.3j 



j=i j=i 

j^k 



n s\n{(Oj - e^-)/2) 
j=i 

J^k 



because 



1 - = -2/sin(/?/2), 



(■ 1 - 1 ) 



and ^'j-^(i^j — &j) = T. The formula for the weights now follows by substituting (4.4) with P := t 
into (4.3). □ 



We can now state the inverse spectral problem and its solution. 



Problem: Given two sets of n mutually interlacing points on the unit circle {Aa:}JJ_| and 
determine the unique unitary llessenberg matrix II = //(71 , . . . , 7„_i . 7^^) and |a| = 1 such that 
the spectrum of II is and the spectrum of //(o7i , . . . , a7„_i , 07^) is 

Solution: Let a := 11^=1 (^^^VA/c)^ calculate the weights by Proposition 4.3, and use the inverse 
unitary QR algorithm to construct II — 7/(7], . . .,7n-i'7n)- 
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